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Abstract We analyze the analytic structure of the Covariant Spectator Theory (CST) contribution to the self- 
energy amplitude for a scalar particle in a theory. To this end we derive dispersion relations in 1+1 
and in 3+1 dimensional Minkowski space. The divergent loop integrals in 3+1 dimensions are regularized 
using dimensional regularization. We find that the CST dispersion relations exhibit, in addition to the usual 
right-hand branch cut, also a left-hand cut. The origin of this "spectator" left-hand cut can be understood in 
the context of scattering for a scalar </) 2 ^ 2 -type theory. If the interaction kernel contains a linear confining 
component, its contribution to the self-energy vanishes exactly. 
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1 Introduction 

The energy spectrum and internal structure of mesons are influenced by both the Dirac structure of the effec- 
tive confinement interaction and the form of the low energy effective one-gluon-exchange interaction, both of 
great theoretical interest. Knowledge of meson spectrum and meson wave functions allows us to unfold in- 
formation about these effective interactions, and is therefore of fundamental importance to hadronic physics. 
Knowing the wave functions also permits the calculation of decay rates and meson transition form factors, 
both of which play an important role in a large variety of reactions, including hadronic contributions to the 
anomalous magnetic moment of the muon [ 1| that is being investigated for possible signs of physics beyond 
the Standard Model, and in the production of lepton pairs whose precise knowledge is mandatory in the search 
for the quark-gluon plasma [2J. This strong interest has lead to the development of current and future experi- 
mental programs to measure meson spectroscopy at Jefferson Lab (CLAS12 and GLUEX programs), CERN 
(COMPASS), and GSI (PANDA at FAIR). Although the principal goal of these experiments is to find "exotic" 
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mesons, such as glueballs and hybrids, there is still much to be learned about mesons that can be described as 
pure qq states. 

Our goal is to construct a unified and consistent relativistic model for mesons as qq bound states that can 
be applied to both heavy and light mesons. The model should satisfy a number of important requirements: (i) 
the formalism used should be covariant, which is essential for the description of systems composed of light 
quarks; (ii) the quarks should be strictly confined; (iii) the structure of the constituent quarks themselves, 
i.e., their self-interaction, should be described consistently through the same confining interaction that acts 
between pairs of quarks; and (iv) the model should reflect the requirements of chiral symmetry (i.e. when the 
bare mass of a quark q\ approaches zero, a massless qtqt bound state must emerge). 

In regard to the covariance property, an important difference between our approach and the very successful 
Bethe-Salpeter Dyson-Schwinger calculations [3 ] is that the latter use a Euclidean metric whereas we stay in 
Minkowski space. Performing a Wick rotation [|4] of the relative-energy component to the imaginary axis, and 
thereby effectively turning the physical Minkowski into a Euclidean metric, has the computational advantage 
that the propagator singularities located on the real axis are avoided. This method thus simplifies the calcu- 
lation of binding energies. However, the Bethe-Salpeter bound-state amplitudes are obtained for unphysical 
energy components. For the calculation of electromagnetic form factors, the Bethe-Salpeter amplitude has to 
be integrated over physical loop momenta, but the reverse Wick rotation crosses singularities of the integrand. 
In addition, for non-zero momentum transfer, the Bethe-Salpeter amplitude has to be known in a frame with 
non-zero total momentum. In principle it can be obtained through a boost from the rest-frame amplitude. 
However, if the rest-frame amplitude is known only for imaginary relative-energy values, the boost parameter 
becomes complex, and performing the boost actually requires knowledge of the amplitude in the whole com- 
plex plane. It is not clear how the required analytic continuation could be performed in a reliable manner 
[6113. These difficulties are often ignored, but Refs. H7] showed in simple model calculations that significant 
differences between form factors calculated in Euclidean and Minkowski metric can be observed. 

A different way to deal with the difficulties of the Bethe-Salpeter equation is to simplify it without leaving 
Minkowski space. A number of so-called quasi-potential approaches were constructed for that purpose. They 
all have in common that the dimension of internal loop integrations is reduced from four to three, which makes 
the resulting equations much easier to manage. We use the Covariant Spectator Theory (CST) [8 9, 10j, which 
has been successfully applied in many different few-body systems (for a short recent review see Ref. (TTj ). 
It was also shown in model calculations with a simple interaction between two scalar particles, for which it 
was possible to solve the Bethe-Salpeter equation including a complete kernel of ladder and crossed ladder 
diagrams in a Feynman-Schwinger representation, that the binding energy obtained with the CST equation is 
actually closer to the exact result than the Bethe-Salpeter equation in ladder approximation (the form that is 
almost always used in applications) III 21 . 

In a consistent theory of mesons, the constituent quark mass should be related to the quark-antiquark 
interaction. In other words, the constituent quark mass should be calculated from a bare mass and the dressing 
generated by the same interaction. For the case of the pion, it was shown analytically that in the chiral limit 
of vanishing bare quark mass, the two-body CST equation yields a massless pion solution, while the quarks 
acquire a non-zero dressed mass through dynamical chiral symmetry breaking. However, in the numerical 
calculations performed so far, the constituent quark mass has been treated as a constant II13II141[T51 or as a 
phenomenological function not related to the kernel 1161 . 

From the structure of the two-body bound-state equation, which will be introduced in Sec. [2] one can see 
that the quark mass, which is a function of the quark four-momentum squared, appears inside a loop integral. 
Thus, we have to be able to calculate it for arbitrary — including negative — values of its argument. But before 
we embark on this somewhat ambitious program, we will study a simpler case, namely that of a scalar particle 
dressed by the exchange of another scalar particle and an effective linear confining interaction. A thorough 
understanding of the properties of the scalar quark mass function should provide useful information also 
about the fermionic quark mass function. In particular, it could yield constraints on the interaction models 
with confinement to be used in the realistic case of quarks as fermions. 

This paper is organized in the following way: after this Introduction, Sec. [2] describes briefly the main 
features of our CST quark model. In Sec. [3] dispersion relations for the self-energy of a scalar particle dressed 
through another scalar particle of smaller mass in a ty 2 % theory are derived in 1+1 dimensions and compared to 
the dispersion relations of a four-point ("bubble") amplitude in the corresponding <\) 2 % 2 theory. These calcu- 
lations are extended to 3+1 dimensions in Sec. [4] In Sec.[5]we consider the self energy for linear confinement 
in 3+1 dimensions, and in Sec.[6]we draw our conclusions. 
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2 Quark-antiquark bound-state equation and quark self-energy 

The development of a covariant model of mesons as bound states of constituent quark-antiquark (qq) pairs in 
the framework of the Covariant Spectator Theory (CST) was initiated in a series of papers [13, 14, 15, 16|. In 
this section we will briefly review the most important features of this model and discuss the aspects we plan to 
improve. We begin the discussion by considering the realistic case of two Dirac particles, and later specialize 
to scalar particles. 

The general structure of the CST equations for the qq bound-state vertex function F is displayed graphi- 
cally in Fig. [TJ The interaction kernel V is an operator in the Dirac space of the two quarks. It can be written 
in the general form 

V = £^ 1 (, ' ) ®^'V( , ' ) , (1) 

i 

where (?f refers to a Dirac matrix of type i in the space of particle 1, and describes the momentum 
dependence of that particular type of interaction (generally confinement plus particle exchanges). Assuming 
that particle 1 is the quark with four-momentum p\ and particle 2 the antiquark with four-momentum p 2 , the 
algebraic form of the two-body equation is 

+ 1 / 2^(2^)3 mkuk 2 )(m + fe^'VQfr, ) . (2) 

We use a notation in which a hat over a particle's momentum explicitly indicates that it is on mass shell, e.g., 
p 2 = m 2 , whereas p t is in general off mass shell. Furthermore, P = p\ + p 2 = p\ + p 2 is the total momentum 
of the meson, m the dressed quark mass, and = (in 2 + k 2 ) l l 2 . The propagator for an off-shell quark with 
four-momentum p is 

1 2 M(p 2 ) + f 

S{p) = m Q -f> + £(p)-ie =Z{p ] M 2 {p 2 )-p 2 -ie' (3) 

where the dressed quark mass function M(p 2 ) and the renormalization factor Z(p 2 ) are related to the quark 
self-energy E(p) through 

Z(p)=A(p 2 )+p'B(p 2 ), (4) 



M(p 2 )=Z(p 2 )[m Q +A(p 2 )] , (5) 
Z{p 2 )=[\-B{p 2 )Y\ (6) 

and mo is the bare quark mass. The dressed mass m is the solution to the equation M(m 2 ) — in = 0. 

It is one of the main features of the CST of a two-body system that in intermediate states the heavier 
particle is placed on its positive-energy mass shell. Since we are particularly interested in obtaining solutions 
for light mesons made of up and down quarks which are treated as equal-mass particles, we have to generalize 
this prescription in a way that restricts either particle on mass shell with equal weight. This leads to the system 
of two coupled equations of Fig. [TJ and Eq. ([2]). The two-channel formulation is also necessary in insure 
consistency in the limit when mo —5- 0, necessary for a study of the chiral limit [13|. For the description of 
heavy-light quark systems, the system of coupled equations can be reduced to a single equation, consisting of 
the first equation of Fig. [TJ with only the first term on the rhs. 

For a realistic model of qq bound states, the kernel must include a confining interaction. Nonrelativistic 
potential models and lattice QCD studies suggest that the strength of this confining interaction increases 
linearly with the distance between the quarks, but there is no general consensus about its Lorentz structure. 
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Fig. 1 The coupled CST equations for the qq bound-state vertex function (triangle). The zigzag line represents the kernel, which 
in this paper is the sum of a confining interaction and the exchange of another particle. A cross on a quark line indicates that the 
particle is on mass shell. 



It is not obvious how a linear potential can be used in momentum space, since its Fourier transform does 
not exist. Nevertheless, Ref. 1131 introduced a limiting procedure which allows an exact solution of bound 
states in a nonrelativistic linear-potential in momentum space, and which can be generalized to the relativistic 
case, although such generalizations are clearly not unique and require additional constraints. According to 
this method, the momentum-space version of the nonrelativistic linear potential V(r) = or can be written in 
the form 

V A (q) (2nf8^(q) j W)j , (7) 



Vl(9) = lim 

£->0 



where q is the momentum transfer, and 

Vh{q) = -W^f- (8) 

The subtraction term proportional to the delta function in |7| is crucial. The Fourier transform of Va alone back 
into r-space diverges as £ — > 0, and the subtraction term exactly cancels the singular part. A straightforward 
covariant generalization is the simple substitution q 2 — > —q 2 , which guarantees that a linear confining potential 
is obtained in the nonrelativistic limit. This has advantages, but also disadvantages, which is why other forms 
have also been explored [ 16 1. Which form is the most suitable to be used in the relativistic kernel is still under 
investigation. 

In past CST calculations, the linear potential was chosen to be purely scalar, and a constant vector interac- 
tion was added to adjust the mass scale. Additional interaction terms can be included if needed, but the exact 
form of the kernel will eventually be obtained from a fit to the meson spectrum. 

However, the choice of allowed interactions in the kernel is constrained. For instance, chiral symmetry 
requires that in the chiral limit, i.e., for vanishing bare quark mass, the pion must become massless. The 
dressed mass of the constituent quarks is generated through the dynamical breaking of chiral symmetry. In 
the CST quark model, this can be implemented in a way that is analogous to the Nambu-Jona-Lasinio (NJL) 
mechanism ffTTIl . where chiral symmetry is dynamically broken through point-like fermion self-interactions. 
In our case, the self-interaction is of a more general kind, because it also includes confinement. 

In fact, when the quark mass function is calculated from the self-energy generated by a certain class 
of relativistic qq kernels, in the chiral limit the qq equations automatically have a zero-mass bound-state 
solution 11131 . In general, this can only occur with relativistic kernels that are invariant under chiral rotations. 
Remarkably, it also holds for scalar confining interactions. Although a scalar term breaks chiral symmetry, 
in the chiral limit it completely decouples from the equations and is therefore an allowed component of the 
kernel. 

As already mentioned in the Introduction, we want to eliminate a lack of consistency in previous calcu- 
lations of the structure of light mesons in the CST, where the dressed quark mass was treated as a constant 
parameter. In this work, we calculate the self energy of scalar quarks which allows us to substitute this constant 
quark mass by the correct mass function. 



5 



P-k 



\ 



Fig. 2 The self-energy diagram. 



A convenient form to represent the mass function is in terms of a dispersion integral which reveals its an- 
alytic structure. Here we have to address a more fundamental question that arises in this context: what kind of 
dispersion relations do the self-energy amplitudes calculated in the CST framework satisfy? This is an inter- 
esting question, because by selecting certain poles in the complex relative energy plane and excluding others, 
the analytic structure of the corresponding amplitudes is surely altered, and it is important to understand the 
consequences of these changes. 

In the following, we will begin by deriving the dispersion relations for the self-energy amplitude for the 
simple case of one scalar particle dressed by another scalar particle with a different mass. Before considering 
the 3+1 dimensional case, it is useful to analyze the problem first in 1+1 dimensions where one can obtain 
analytic results for both the Feynman and CST amplitudes without any necessity to regularize divergent loop 
integrals. 



3 Two scalar particles in 1+1 dimensions 

The aim of this section is to demonstrate how the self-energy amplitude is derived within the CST framework. 
We can compare our results with another approximation used by Becher and Leutwyler in Ref. H 1 81 in the 
context of heavy baryon chiral perturbation theory, and we apply some of the methods to derive dispersion 
relations described there. 

The self energy corresponds to the one-loop graph shown in Fig. [2] At first, in order to keep the calculations 
simple and to avoid singularities emerging from divergent integrals, we consider the scalar loop diagram 
involving two scalar particles of masses M and m with M > m in two-dimensional Minkowski spacetime^] 
This simple example serves as a review of the basic concepts and ideas. Then, in the next section we turn to 
more intricate case of four space-time dimensions where the scalar-loop integral diverges, which we regularize 
through dimensional regularization. 



3.1 Dispersion relation 

This subsection reviews the dispersion relation for the self-energy diagram in two dimensions and serves as a 
preparation for the subsequent sections. From Fig.|2]the familiar Feynman scalar-loop integral is given by 

'"'"'/Wis 

with 

A x =M 2 -k 2 -i£x, 

A 2 =m 2 -(P-k) 2 -ie 2 (10) 

where k^ = (k°,k l ) is the internal two-momentum and P^ = (P°,P l ) is the total two-momentum with Man- 
delstam s = P 2 . The values of coupling constants are irrelevant in this context and thus are set to 1 throughout 
this work. The 'ie' prescriptions in the denominators of the propagators impose the boundary condition that 
positive-energy states propagate forward and negative-energy states backward in time. Note that we have 
equipped the two masses with different imaginary parts whose zero limits are taken one after the other. Al- 
though this requirement is not needed for now it will become relevant in the next subsection. 



We use the standard convention for the metric tensor diag(g fiv ) = (1,-1). 
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Fig. 3 The real (solid) and imaginary (dashed) parts of E(s) for the choice M = 2 and m = 1. 



The product of the two propagators in Eq. |9) can be rewritten using the standard Feynman parametrization 



as 



1 



dz 



A X A 2 Jo {A 2 z+A l (l-z)) 2 
Performing the integrations over d 2 k gives 

1 f 1 dz 



E(s) = 



with 



C(z,s)=z(l-z) 



An Jo C{z,s) 
M 2 



1 



£2 
1-Z 



(11) 



(12) 



(13) 



In the z-interval of integration [0,1] the imaginary part in the square brackets is always negative, therefore 
it is safe to make the following replacement in that region without changing the 'ie' prescription: 



£2 | £i 
1— z z 



£ . 



The integral of Eq. ( 12 1 is elementary and can be written as 



= - 2n p(s + ie) ^ arCtan I- 11- ^ + + arCtan ^ + ^ + 



where 



p(s) = y/(s-S+)(s-S-) 

is the familiar two-body phase space factor with 



s + = (M + m) , S- = (M — m) and T}±(s) = — i 



. s±(M 2 -m 2 ) 



(14) 

(15) 
(16) 

(17) 



The real and imaginary parts of E(s) are depicted in Fig. [3] 

Note that the imaginary part of E(s) is non-zero only above the threshold 5 > s + where it is kinematically 
possible that both particles are on-mass shell. It has a singularity at 5 = s + and a discontinuity along the branch 
cut on the positive real axis from s + to °° as illustrated in Fig.|4j The 'ie' prescription in the arguments ensures 
that ImZ(i) is always evaluated on the correct side of the cut, which is in the case at hand the upper-half 
complex s plane. We proceed with our discussion by showing that the self energy E(s) satisfies a dispersion 
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relation. In order to find the corresponding integral representation we perform a variable transformation of the 
z integral in Eq. ([T2| defined by the map 



it \ 

s(z) = y 

The inverse transformation involves two roots: 



M 2 

z 



(18) 



z ± (s') = — [s'+M 2 -m 2 ±p(s')] . 



(19) 



On the interval of integration < z < 1, the negative branch covers the interval < z_ < z(s+) =Mj (M + m) 
and the positive branch covers the interval z(s + ) < z+ < 1. Furthermore, each of these intervals is ma pped 
onto the same s' interval s + < s' < °° and this interval is thus covered twice by the mapping of Eq. ( 18 1. 
Written in terms of the new variable s', the integral of Eq. ( 12 1 acquires the expected dispersive form 



1 

2n 



ds' 



1 



p(i') s' — s — ie 



it 



'ds'lmZ{s') 



ie 



(20) 



with 



ImEM 



2p{s) 



-eis 



(21) 



It can be checked numerically that the analytic expression of Eq. ( 15 I agrees with the dispersion integral 
of Eq. pi. 



3.2 Dispersion relation in the spectator theory 

In this section we turn to the CST framework and derive the corresponding dispersion relation for the self 
energy. Applying the CST program to the self-energy diagram means that the heavier particle of mass M is 
placed on its positive-energy mass shell. An elegant procedure to work out this 'spectator' contribution has 
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been presented in Refs. [9 18 , 19 1. The idea is to split up the product of the two propagators of Eq. (Ill into 
two terms according to the algebraic identity 



1 1 1 



(22) 



A X A 2 Ai(A 2 -Ai) A 2 {A 2 -A X )' 
or, equivalently, in terms of the z integration as 

f 1 =r r * m) 

Jo (A^+Aiil-z)) 2 Jo (A 2 z+A,{\-z)) 2 J i (A 2 z+A,{\-z)f K ' 

We can design this representation in such a way that the desired spectator contribution is given entirely and 
only by the first term on the right-hand side of this equation by choosing an appropriate "e-prescription". In 
order to see this explicitly we go to the overall rest frame, where P = (W, 0) with W = y/s. Inserting for the 
propagators this contribution becomes 

Csn ; J (2k) 2 (E k -ko-i£i)(E k + k -i£i)[m 2 -M 2 -s + 2Wk -i(e 2 -e 1 )] ' v ; 



where E k = yM 2 +k 2 . At this point it is obvious why we have equipped the two masses with different 

imaginary parts: it simply avoids the cancellation of the two imaginary parts and thus provides a well-defined 
'ie' prescription of the pole in the square brackets. Otherwise, the absence of an 'ie' prescription would make 
the integral over k ambiguous [ 1 8 1 . In order to obtain the spectator contribution of the heavy particle on its 
positive-energy mass shell we assume the following 'e-limiting' prescription: first take the limit £i — >• and 



then take the limit £ 2 — > (such that the difference £ 2 —E\ in ( 24 1 is always positive) and assume that W > 0. 
In this case there is only one pole in the lower-half complex ko plane, the spectator pole at ko= E k — \£\ . Now 
we can carry out the ko integration and we recognize the familiar CST term ( see, e.g., Ref. ffTTIl ): 

£CST(S) = - / ^ F 2 1= — TyT (25) 

./ 271 2E k [e{-{W -E k Y\ 



where e k = ym 2 +k 2 . 

Note that we would have obtained the same result if we would have taken the opposite 'e-limiting' pre- 
scription provided W < 0. For the other two possible assumptions the single pole would have been in the 
upper half plane at ko = —E k + i£i , corresponding to taking the negative-energy pole of the heavy particle. 
Therefore, the 'e-limiting' prescription together with the sign of W simply fixes the position of the pole 
ko = (s +M 2 — m 2 )/2W in the complex plane. 

With the same 'e-limiting' prescription as above and in an analogous manner one can show that the 
second term on the right hand side of Eq. ( |22"| ) corresponds to the negative-energy pole contribution of the 
light particle of mass in. 

In order to find an analytic expression for IcstO 5 ) we consider the first term on the right hand side of 



Eq. (23 I and integrate over d z k to obtain 



W*) = ~/"^r- (26) 

47T./0 C(z,s) 

This resembles the result in Eq. ( |12[ ), with the difference that the upper limit is now °o instead of 1. Before 
performing the z integration we have to analyze the 'ie' prescriptions over the relevant s regions in order to 
obtain the correct sign of ImZcsT(^)- With our 'e-limiting' prescription defined above the imaginary part in 



the square brackets of Eq. ( 13 I is negative for < z < 1 and positive for z > 1. This, in effect, results in a 
sign change of the 'ie' prescription on how to treat the pole at s'(z) = s. Note that in the interval < z < 1 
the function s'(z) has a minimum s + at z = M/ (M + m) and approaches °° as z — > + or z — > 1 ~. Further, in 
the interval 1 < z < °° the function s'(z) has a maximum i_ at z = Mj [M — in), approaches — °° asi-> 1 + 
and becomes as x — > °°. Therefore we have a '+ie' prescription in the interval s+ < s < °° and a '— ie' 



prescription in the interval — °° < s < S-. We can now perform the z integration in Eq. (26 1 to obtain 

-S^heylf +arctan[77+( S -ie)]}, -o° <s<s 



-CST 



[') = < 



2^(1 +arctan[T] + ( 5 )]} , S-<s<s+ (27) 

2Hp^{l + arctan[77+( S + ie)]}, s+ < s < o. . 
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Fig. 5 The branch cut discontinuity of ImZcsT ( J + i£) for the choice M = 2 and m=\. 



Our next task is to write Zcst^) in a dispersive representation similar to Eq. (20 1. To this end we proceed 
in an analogous manner as in Sec. 3.1 by writing Eq. (26 1 in terms of the s'(z) variable. Note that the z 
interval 1 < z < z(s~) covered by the negative branch Z-{s) is mapped onto the s' interval — °° < s' < S- and 
the interval z{s-) < z < °° covered by the positive branch z+(s) is mapped onto S- > s' > 0. Therefore, the 
interval < s' < s_ is covered twice by the map z — > s'(z). As the analysis after Eq. (26 1 revealed there is a 
'— ie' prescription in the interval — °° < s < S- and a '+ie' prescription in the interval s + < s < °°. Taking this 
into account we obtain 



-CST 



1 

2tz 



1 /■» 

2 J 



Jo 



d.v' 



1 



p{s') s' -s+ie 



r ds' l 1 

J s+ p-(?)7^iej (28) 



from which we find the following expression for IiilEcstM 



ImXcsT(^) = 



1 



2p(s) 



~0{-s) + 0(s)0(s- -s) + 6(s-s+) 



(29) 



Note that ImZcsT(^) vanishes in the interval S- < s < s + , it has a discontinuity at s = and it is singular at 
s = i_ and s = s + . Further, it has two cuts along the real axis, the usual 'right hand cut' going from s + to °° 
and a 'left-hand cut' going from — °° to s_ as shown in Fig. [5] The numerical computation of the dispersion 
integral of Eq. (28 1 agrees with the analytic expressions of Eq. (|27[). 



A comparison of Eq. (28 1 with Eqs. (22 1 and (20 1 reveals where the 'spectator left-hand cut' stems from: 
it corresponds to (minus) the negative-energy pole contribution of the light particle that is omitted in the 
spectator approximation. 

We shall now discuss the difference between the spectator approach presented so far and the alternative 
approach given by Becher and Leutwyler in Ref. [18]. Instead of keeping the positive-energy pole of the 
heavy particle, Becher and Leutwyler have rather kept the positive-energy pole of the light particle. The 
corresponding contribution is most easily obtained by simply interchanging A\ with A2 on the right hand 
sides of Eqs. ( 22 1 and ( 23 1, which corresponds to interchanging the masses m and M and the momenta k and 
P — k. The relevant term is then given by 



dv 



{A iy +A 2 {\-y)f 



dz 



(A 2Z +A 1 (l-z)) 2 



(30) 
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Fig. 6 The branch cut discontinuity of IitiZbiX 5 + i £ ) f° r the choice M = 2 and m = 1. 



where we have transformed the integral from a y integration to a z = 1 — y integration. This term leads to the 
desired contribution provided the limit £2 — > is taken before the limit £\ — > and W > 0. This alternative 
'e-limiting' prescription contrasts the one from the CST contribution and ensures that there is only one pole in 
the lower ko plane, namely the positive-energy pole of the light particle. With this prescription, the imaginary 
part in the square brackets of Eq. ( 13 1 changes the sign in the interval — °° < z < 1 at z = 0. Further, the 
interval — °° < z < covered by the positive branch x+(s') is mapped onto the interval > s'(z) > — 00 . This 
effectively results in a '+ie' prescription in the interval s + < s < °° and a '— ie' prescription in the interval 
— °° < s < 0. Integrating Eq. ( 30 1 over d 2 k and then performing the z integration yields the analytic result 



£bi» = 



1 

4n 



dz 



• C(z,s) 



= < 



-2Spi=sy{f +arctan[77_( S -ie)]}, -co< s< 
-2^y{f +arctan[T]_(s)]} , < s < s+ 

{| +arctan[?7_(5 + ie)]} , s+ < s < 00 



(31) 



27rp(s+i£ 



The cut structure of IitlEbl^) is displayed in Fig. [6] 

As expected, by interchanging the masses M and m the analytic results for Zcst^) and ■Ebl('S') are mapped 
onto each other, as can be seen explicitly by comparing the expressions in ( |27| l and ( |31| l. Finally, we obtain 
the following dispersion relation for the self energy in the Becher-Leutwyler approach: 



1 

2n 



1 r° ds' 1 

2 7-oo p(s') s'-s + ie 



ds' 1 

J s+ p(s') s' - s - ie 



(32) 



with 



IitlEbl (s) 



2p(s) 



■e(s- 



(33) 



which agrees numerically with the analytic result of Eq. ( |31| . 

The real and imaginary parts of Zcst(^) an d ^bl('S) are compared with the ones of L(s) in Fig. [7] Both 
Ucst(^) an d £bl(s) have left-hand cuts not present in the full amplitude. It is instructive to look at a related 
(j) 2 X 2 theory to better understand the origin of these cuts. 
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(a) 



(b) 



Fig. 7 The real parts (a) and imaginary parts (b) of £cst(- s ) (dashed line) and ZrlM (dotted line) are compared with L(s) (solid 
line) for the choice M = 2 and m=\. Note that in (b) for s < ImZcsT(^) and ImXBL(s) lie on top of each other and for s > s + 
all three lines lie on top of each other. 



3.3 The 'left hand cut' in the two-particle scattering amplitude 

From the topology of the self-energy graph in Fig. [2] with one incoming and one outgoing line it is obvious 
that there exists only one channel, the s-channel we have discussed so far. If one considers, however, elastic 
scattering of two particles, then there are two topologically different Feynman diagrams that contribute to the 
second order amplitude. These are the s- and the u-channel diagrams which are related by crossing symmetry. 
In this section we will show that for two-body scattering it is the u-channel contribution that gives rise to a 
'left hand cut' along the real s-axis similar to the spectator 'left hand cut' discussed in the previous section. In 
order to see this explicitly we consider a </) 2 ^ 2 -type theory in 1+1 dimensional Minkowski spacetime of two 
interacting scalar particles of masses M and m. The corresponding second-order contributions to the scattering 
amplitude are the s- and the u-channel bubbles shown in Fig. [8j 



P-p' 



P-p. 



P-k 



P-p'-p + k\ \k 



P-p 



P 



(a) 

Fig. 8 The s-channel (a) and the u-channel (b) bubble diagram. 



(b) 



The s-channel bubble contribution B{s) of Fig.|8]coincides with E (s) of Eq. ^ for the self energy which 
has already been discussed in Sec. |3.1| with the analytic result and the dispersion relation given by Eqs. ( |15| l 
and (20 1, respectively. Therefore, we concentrate in this section on the u-channel bubble B[u) of Fig.[8](b). It 
is given by 



f d 2 k 1 



(34) 



with 



B 2 =m 2 -{P-p' -p + kf -ie 2 



(35) 
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and with Mandelstam u = (P — p' — p) 2 . By introducing a Feynman parameter and performing the d 2 k inte 



gration, this contribution can be put into a form that is equivalent to the one for L(s), see Eq. ( 12 1: 



1 f 1 dz 

B ( u ) = -J w rY^V (36) 
4n Jo C(z, u) 

The fact that B{u) and B{s) have the same functional form is a consequence of crossing symmetry. Conse- 
quently, B(u) and B(s) satisfy the same dispersion relations: 



„, . 1 f°° ds'lmB(s') 
% J s+ s'-u—ie 



(37) 



with 



]mB(u) = -- T ^G{u-s + ). (38) 
2p{u) 

If we consider only cases where Mandelstam f is small, we can use the relation 

s + u^lM 2 +2m 2 (39) 

to express u in terms of s, which gives us an approximate expression for the u-channel contribution as a 
function of Mandelstam s: 

B(u) « B(2M 2 +2m 2 - s) = B u (s) . (40) 



With this approximation the dispersion relation for B u (s) obtained from Eq. (37 1 reads 



1 f s - ds"lmBJs") „ , 



B u (s) = — —r, ^V^ with s" =2M L +2m z - s (41) 

and 

JmB u (s) = - T ^ r -B(s--s). (42) 
We see that the u-channel bubble contribution expressed in terms of the Mandelstam 5 variable has a 'left 



hand cut' along the real s-axis going from — °° to i_ . With the help of Eq. ( 39 1 we can approximate the total 
second-order contribution as a function that depends solely on s: 

B(s) +B(u) w B(s) +B u {s) = J£ % {s) . (43) 

The CST and Becher-Leutwyler contributions to the s-channel bubble of Fig. [8j Scst(*) and Bbl(s), are 
simply the same as for the self-energy diagram of the preceding section, Zcst(s) an d £bl(s), respectively. 

In Fig. [9] the real parts of j%%{s), Bcst{s), Bbl(s), B(s) and B u (s) are shown for different mass ratios 
a = M/m. First, note the symmetry of the total contribution ^2{s) about the point s = s m = M 2 +m 2 as con- 
sequence of crossing symmetry and l«0. Second, from Fig. [9] we observe that for the case of equal masses, 
i.e. a ~ 1, #cst(s) anc ^ ^bl(s) become identical, as can be seen explicitly from the analytic expressions ( |27| > 
and ( |31| l, which differ only by interchanging the masses M <R- m. Third, Fig. |9]reveals that for the case when 
M becomes very large compared to m, i.e. a 3> 1, the Becher-Leutwyler contribution Sbl(^) goes over to the 
s-channel contribution B(s) and the spectator contribution Scst(^) approaches the total contribution j^%{s). 
Actually, the latter observation is simply the confirmation of the CST cancellation theorem for </> 4 -type theo- 
ries discussed in Ref. 1111 : in the limit a — > °° the sum of the s- and u-channel diagrams becomes identical to 
the part of the s-channel diagram where the heavy particle is placed on its positive-energy mass shell. There- 
fore, we conclude that for mass ratios a S> 1, the Becher-Leutwyler contribution is a better approximation 
to the self-energy diagram (or equivalently the s-channel bubble), whereas the CST contribution is a better 
effective description of the total bubble diagram (s- plus u-channel bubble). Consequently, for K>1 the 
difference between the Becher-Leutwyler result and the CST result for the bubble diagram is the u-channel 
bubble contribution. 
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Fig. 9 The real parts Re^O) (solid), ReBcsT (*) (dashed), ReB B LM (dotted), ReB(s) (long-dashed) and ReB„ (s) 

(dot-dashed) are shown for m = 1 and for different choices M = 1.05 (a), M = 1.5 (b), M = 2 (c) and M = 9 (d). 



4 Two scalar particles in 3+1 dimensions 

We now apply the above ideas to the physical 3+1 dimensional Minkowski spacetime. In this case we are 
confronted with divergent integrals for the self-energy and for the bubble-diagram contributions. In order to 
find the dispersion relations for such diagrams we have to separate the regular from the singular parts. This 
will be achieved by a procedure [ 18 ] based on dimensional regularization. 



4.1 Dimensional regularization and dispersion relation 

Let us first analyze the scalar-loop integral for the self energy of Fig. [2] in d spacetime dimensions (with d 
being an arbitrary value) and then perform the limit d — > 4: 



d d k 1 

(2n) d AiA 2 



(44) 



Here all two- vectors in A\ and A 2 of Eq. ( 10 1 have been replaced by d- vectors]^] Introducing one Feynman 
parameter and performing the d d k integrations using the standard techniques leads to 



Lis) =-r 2-- 



dz 



2 J Jo ( 4jc )S 



(45) 



2 As in Sec[5] we use the standard convention for the metric tensor diag(g A , v ) = (1, — 1,..., — 1) with jtt, V = 0, . . . ,d— 1. 
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This expression is singular for d = 4 due to the F-function in front of the integral which reflects the logarithmic 
divergence of the loop integral (44 1 in four dimensions. In order to separate the finite from the singular parts 



of E(s) we write d = 4 — <5 for 5 > 0. Then, by expanding [C(z,s)] 1 in \45\ about small values of 8 we 
obtain 



2F 1 



8{4n) 2 



1 



dzlnC(z,s) 



(46) 



Here the first term in the square brackets is identified as the singular part as 8 — > 0. Since this part is indepen- 
dent of 5 we define the finite part of E(s) in d = 4 dimensions by a single subtraction at a point s = sq: 



dzln 



C(z,s ) 



This integral is elementary and can be expressed as 



E(s) 



1 fipO + ie) 



{Alt) 



{arctan [tj + (j + ie)] + arctan [tj_ (5 +ie)]} 



ippo + ie) 



{arctan [tj+ (j + ie)] + arctan [tj_ (j + ie)]} 



In 



M 



[zo(j+ie)-zo(jo+ie) 



(47) 



(48) 



with 



z (s) = -(s+M 2 -, 
2s 



(49) 



Note that the '+ie' prescription again ensures that the imaginary part ImZ(i) is always evaluated on the 
correct side of the branch cut where it has a discontinuity, which is in this case always on the upper-half 
complex 5 plane. The subtracted dispersion relation for E(s) follows from Eq. \4l\ by an integration by parts 



and by transforming the integral from z to s'(z) (see the discussion in Sec. 3.1 

s — So /'"" ds'Im.£(.s') 



Z(s)=Z(s )- 



K 



•ie 



with 



hn z(s) = -^-e(s-s + ). 

Ions 



(50) 



(51) 



This dispersion relation resembles the corresponding one in two dimensions, see Eq. (20 1, with a branch- 
cut discontinuity along the real 5 axis from s + to °°. 



4.2 Spectator dispersion relation 

The spectator contribution Zcst(^) to the self-energy diagram in d dimensions is obtained from E(s) in the 
same manner as in Sec. |3.2| by splitting up the product of the propagators according to Eq. (23 I. This gives 



tcsr(.v>: -r[2-i) r^d- 2 (z,s). (52) 

(4tt)2 



2 



Similar to £(s) of the previous section this contribution is singular in d = 4 dimensions because of the F- 
function in front of the integral. Furthermore, for d = 4, the z integral does not converge anymore due to the 
upper limit at z = °°, which makes the CST case somewhat more intricate as compared to E (s) from before. 
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In order to work out the finite part of Zcst{s) we follow the procedure presented in Ref. 0181 . To this end we 
note that C(z,s) is a second-order polynomial in z that can be written in the form 



C(z,s)=C {s)+s[z-zo(s)} 2 
where (suppressing the small imaginary parts) 



Co (*) = 



p 2 M 

4s 



Then 



Ci- 2 (z,s) = C (s)C-2- 3 (z,s) +s [z -zo(s)} 2 Ci- 3 (z,s) 



(53) 



(54) 



(55) 



where the second term is proportional to the derivative of C2 2 (z,s) with respect to z. An integration by parts 
gives 



-CST 



(s) = 



1 



2T{3- 



(4tt)i (d-3)(d-4) 



[z-zo(s)}d- 2 ( z ,s) 



+ (d-4)C Q (s) rdzC2-\z,s)\ . (56) 
=0 Jo J 



The first terms in the curly brackets are singular for d = 4. Let us concentrate on the remaining finite term of 
the integral. In four dimensions this integral converges and gives the finite part of the CST contribution (for 



-CST 



2C (s) 



dz 



(4tt) 2 Jo C(z,s) 



'-^{f + arctan^-ie)]}, -°°< s<s - 
"(^{f +arctan[T] + (i)]} , (57) 
- i ^{f + arctan[T J+ (,+i £ )]} ) s + <s<~ . 



Note that this function has a pole at s = 0. Following [18], this suggests the twice- subtracted dispersion 
relation 



-CST 



, z (s) - p 2 
zo{s+) 47 



ri r° 




2I.J 


-f 

Jo 



ds' 



1 



p(s') s' -s + ie 



[°° ds' 1 1 

X + pl7)7^ii/' (58) 



from which the imaginary part follows as 



ImEcsT(^) 



16ns 



d(-s) + e(s)d(s--s) + 6(s-s+) 



(59) 



Thus, we encounter the typical CST branch-cut structure of the two-dimensional case, see Eq. ( |28| l. 

The Becher-Leutwyler contribution £bl( s ) to the self-energy diagram in four dimensions is expressed by 
similar formulae (see Sec. 3.2 1 and will not be given here explicitly. 



4.3 Bubble diagram 



The contributions to the s-channel bubble diagram (Fig.[8](a)) in four dimensions, B(s), Bcst(s) and 5bl(s) 
coincide with the ones for the self energy, E(s), Zcst( s ) an d £bl(s), respectively. What is still missing for the 
calculation of the full bubble diagram in four dimensions is the contribution of the w-channel bubble diagram 
of Fig. [8] (b). In the approximation of Eq. (40i this contribution can be written as a function of Mandelstam 
5 only. Using the methods of Sec. 4.1 this contribution satisfies the following once-subtracted dispersion 
relation: 



B u (s) =B u {s Q )+B u (s) 



-so r 



ds'lmB u (s') 

(s'-s )(s'- 



(60) 
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Fig. 10 The imaginary parts Irrk^2 ) (solid), ImBcsT(s) (dashed) and IitiBbl(-s) (dotted) are shown for the choice M = 2 and 
m = 1. Note that for s > s_ all three curves, and for s < ImBcsTM and IitiBblCs). lie on top of each other. 



with 

\mEJs) = - — r0(*_-*). (61) 

v ; 167r(2M 2 + 2m 2 -^) v ' 

As expected from the result in two dimensions, we see that also in four dimensions the M-channel bubble 
contribution gives raise to a 'left-hand branch cut' in the complex s-plane along the real axis from s = — °° 
to s = S-. Finally, we can write down the approximate expression for the finite part of the total second-order 
bubble amplitude: 

Z# 2 (s)=B(s)+B u (s). (62) 

The imaginary parts of the contributions are compared in Fig. [10] 

It is interesting to see the difference between the imaginary parts in two and in four dimensions. While 
the former are singular at 5 = s± (due to the p function in the denominator), the latter, in contrast, vanish at 
s = s± (p function in numerator). Furthermore, in four dimensions the imaginary parts of BcstC 5 ) and 5bl(s) 
are singular at s = and become positive for s <0, whereas in two dimensions they are always negative. 
In order to compare the real parts we choose the subtraction point for B(s) and B„(s) to be at threshold, 

for different mass ratios a = tr ■ As in two 
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i.e. sq = s+. The real parts of the contributions are shown in Fig. 
dimensions, the symmetry of ^2 {$) about the point s m = M 2 + m 2 is a consequence of crossing symmetry. 
Further, as expected from the symmetric appearance of the masses in the expressions, for equal masses (a ~ 1) 
ReficsT^) an d ReBBiX?) become identical. Also, for a — > °° the CST result approaches the total bubble 
contribution. This happens, however, not as fast as in two dimensions. Nevertheless, we can still conclude 
that the CST result approximates the sum of the two bubble diagrams quite well for M S> m. Interestingly, in 
four dimensions the Becher-Leutwyler result does not converge to the total self-energy contribution as a — > °° 
and is therefore not as good an approximation in this limit as it is in two dimensions. 



5 Confining interaction between scalar particles 

We consider now a simplified model of mesons as bound states of scalar particles in the framework of CST, 
whose interaction consists of a confining part and another one that is described through the exchange of a sec- 
ond scalar particle. In order to make this description consistent, the self-energy of these scalar particles, which 
we will call simply "scalar quarks", has to be calculated from the same interaction kernel. The dressing by 
the scalar exchange has already been calculated in Secs.[3]and|4] so we only need to calculate the contribution 
from the confining part. 
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Following the general form given in Eqs. and (jHJl, a relativistic generalization in CST can be 

achieved by H 1 31 



VRL{P,k) = lim 

£->0 



V^pft-lEkilnfd^Xp-k) 



3/,/ 



d 3 £ 



2E k ,(2n) 



and 



with 



Vra(M) 



87TCT 



(-^ 2 + e 2 ) 2 ' 



The expression for the self-energy of a scalar quark dressed by Vrl in 3+1 dimensions is then simply 

d 3 k 



'L.CST 



2^(2?r) : 



;Vrl(M)=0, 



(63) 

(64) 
(65) 

(66) 



which is easy to verify from the form of the kernel ( |63[ ). That the integral of Vrl should vanish follows from 
the fact that Vl(<?) of Eq. |7} is the Fourier transform of the linear potential V (r) = or, and the inverse Fourier 
transform must yield zero at r = 0. This yields the condition 



d 3 q 



Vl(«)=0, 



(67) 
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and Eq. (66 1 can then be seen as a covariant generalization of (67 1. 

To summarize, we have shown that the linear confining interaction does not contribute to the scalar quark 
self-energy at all. If the total kernel consists of a confining part and another scalar exchange as described 
in Sec. |4j then the total result for the self-energy is given by Eq. (58 1. While the scalar confinement makes 
no contribution to the self energy, it does contribute to the structure of the meson. This can be understood 
physically if we think of the self energy as a correction that occurs at short distances (i.e. at r = 0), while the 
meson wave function is sensitive to the kernel at all values of r. 

For completeness it should be mentioned that in 1+1 dimensions the same result is obtained. 



6 Summary and Conclusions 

In this paper, we continue the development of a unified and consistent covariant quark model for mesons in the 
framework of the CST that was initiated by Gross and Milana [13 14 15|. The main aspect of this model we 
aim to improve is the consistent implementation of the quark mass function, in the sense that it is calculated 
from the same interaction that binds the two quarks together. In a first step, we investigated the simpler case 
of scalar quarks and calculated their self-energy amplitude, obtaining a representation for the mass function in 
form of a dispersion integral. The CST approach modifies the analytic structure of field-theoretic amplitudes 
by selecting certain pole contributions and neglecting others. It is therefore interesting to see that the CST 
amplitudes actually do satisfy dispersion relations. 

To avoid complications due to divergences in loop integrals, we performed these calculations first in 
1+1 dimensions. We found that the CST amplitude satisfies a dispersion relation which has a left-hand cut, in 
addition to the usual right-hand cut. This is similar to an approximation investigated by Becher and Leutwyler, 
which also has a left-hand cut. In order to understand the origin of the left-hand cut in the CST amplitude, 
we calculated also the dispersion relations for the scattering amplitude of two scalar particles in a <\) 2 % 2 - 
type theory to second order. The calculation in the s-channel is essentially the same since only the external 
particle momentum has to be written as the sum of two particle momenta, and therefore the same cut structure 
is obtained. However, there is also a u-channel which generates a left-hand cut. Thus, the total amplitude 
has both a left-hand and a right-hand cut, just as the s-channel CST amplitude with the heavier particle on its 
positive-energy mass shell. This is a natural result, since one of the basic ideas of the CST is that its amplitudes 
approximate the sum of direct and crossed diagrams, which in this case means s- and u-channels. On the other 
hand, in the case of the self-energy amplitude at one-loop level, there is no u-channel, and therefore a simple 
physical interpretation for the presence of a left-hand cut is not obvious. At this time the interpretation and 
treatment of this left hand cut is an open question to be discussed in future work. Repeating these calculations 
in 3+1 dimensions using dimensional regularization gave the same cut structure for the self energy. 

Because of the CST constraint on the confining interaction, the self-energy of a scalar quark due to a 
linear confining interaction vanishes exactly, and does not contribute to the mass function at all. The confining 
interaction will contribute to the structure of the meson because the meson wave function is sensitive to the 
kernel at all values of r (or q) while the self energy is a correction that occurs only at short distances where a 
purely linear confining interaction is zero. It would be amusing to study a simple toy model of scalar particles 
bound solely through a linear confining interaction where a constant mass for the constituent particles would 
be part of an exact and consistent solution of the system. We postpone discussion of this toy model for another 
time. 

We emphasize that the vanishing of the self energy contribution due to confinement holds only for scalar 
particles and does not carry over to spinor quarks. With spinor quarks the structure of the confining interaction 
can be divided into a part that gives no contribution to the self energy but may violate chiral symmetry (and 
therefor decouples in the chiral limit) and a chirally invariant part that will, in general, contribute to the quark 
self energy. These issues will be discussed when we extend these studies to the case of spin- 1/2 quarks with 
more realistic interactions, including confining interactions with a spin dependent structure, and when we 
consistently calculate the mass function, the meson spectrum, and the meson wave functions. 
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